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ABSTRACT 

Quasars are powered by accretion onto supermassive black holes, but the problem of the duty cycle related to 
the episodic activity of the black holes remains open as one of the major questions of cosmological evolution of 
quasars. In this Letter, we obtain quasar duty cycles based on analyses of a large sample composed of 10,979 
quasars with redshifts z < 2.1 from the Sloan Digital Sky Survey (SDSS) Data Release Three. We estimate 
masses of quasar black holes and obtain their mass function (MF) of the present sample. We then get the duty 
cycle S(z) = 10"^ ^ 1 based on the Soltan's argument, implying that black holes are undergoing multiple episodic 
activity. We find that the duty cycle has a strong evolution. By comparison, we show that evolution of the duty 
cycle follows the history of cosmic star formation rate (SFR) density in the Universe, providing intriguing evidence 
for a natural connection between star formation and triggering of black hole activity. Feedback on star formation 
from black hole activity is briefly discussed. 

Subject headings: black hole physics — galaxies: active — galaxies: evolution — galaxies: nuclei — quasar: 
general 



1. INTRODUCTION 

Supermassive black holes are relics of quasars in the Uni- 
verse (Soltan 1982; Rees 1984, 1990). Evolution of quasars 
is led by switching on and off accretion onto the black holes. 
During their entire evolution, how many times and how many 
black holes are triggered? what is the triggering mechanism? 
and why do quasars switch off? 

The duty cycle, defined as the fraction of active black holes 
to their total number is a key to tackling the above problems 
(Richstone et al. 1998; Martini 2004). A popular method to get 
the duty cycle invokes the continuity equation and the MF of 
quasar black holes. With an assumption that quasar black holes 
have the same Eddington ratio, their MF can be obtained from 
the luminosity function and finally the duty cycle can be found 
from the continuity equation (Small & Blandford 1990; Mar- 
coni et al. 2004). This is a convenient way to discuss evolution 
of the black holes, but the degeneracy of the Eddington ratio 
and the duty cycle still holds. Actually, not only are the Edding- 
ton ratios not constant for different quasars at different epochs, 
' but they appear to be quite scattered (Vestergaard 2004). The 
duty cycle is poorly understood as a statistical parameter tracing 
the evolution of quasar populations. 

In recent years, there has been much progress in estimating 
black hole masses both in nearby galaxies and distant quasars. 
The empirical relation of reverberation mapping allows us to 
conveniently obtain the black hole masses from a large sample 
and directly get their MF. Thus it becomes realistic to get new 
clues to understand the evolution of quasars. Fortunately, by 
invoking the MF, we can decouple the degeneracy of the duty 
cycle and the Eddington ratio to get the duty cycle. 

In this Letter, we use available SDSS data to directly get the 
MF of the black holes so as to discuss the duty cycle based on 
the MF and find that there is a strong cosmological evolution. 
Our calculations assume a cosmology with the Hubble constant 
Ho = 70 Mpc"' km s"', = 0.3, and = 0.7. 

2. BLACK HOLE EVOLUTION 
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Active black hole evolution can be described by the MF 
<i>(M,,z), which is defined as $(M,,z) = d^N/dM,dV, where 
M, is the black hole mass, d^V is the number of quasars within 
the comoving volume element dV and mass interval dM,. The 
number density of quasar black holes is then given by Nq^oiz) = 
$(M,,z)dM, at redshift z, where M* is the flux-Umited 
black hole mass in a survey. If N{Mt^z) is the MF of all black 
holes including the active and the inactive ones, the total num- 
ber density of black holes is A^tot(z) = Xw. M{M,,z)AM, at red- 
shift z- The relation between N{M,,z) and <I>(M,,z) can be 
written as (Small & Blandford 1992; Marconi et al. 2004), 

$(M.,z) = ^(M.,z)AA(M.,z). (1) 

where the duty cycle, 5{M,,z) is a function of black hole 
mass and redshift. Integrating equation (1) over all black hole 
masses, we have a mean duty cycle in terms of their number 
density 

/^.$(M.,z)dM. _A^^^„(^) 



Sxiz)-- 



J^,AfiM„z)dM, Mo,(z)' 



(2) 



This duty cycle represents the relative number of the active 
black holes to the total. Multiplying by M, and integrating Eq. 
(1) over all black hole masses, we have an averaged duty cycle 
weighted by the masses of the black holes 

/^.$(M.,z)M.dM. (M^iz^ 

^'^'^ = j,,M(M„z)M,dM, = -wm:;^^'^'^' 

where (M,(z))q,o = «'(M.,z)M.dM./A^qso is the averaged 
mass of the quasar black holes at redshift z and (M,(z))aii = 
A/'(M,,z)M,dM,/A^tot for all of active and inactive black 
holes. We show below {M,(z))qso = {M,(z))m forM, >Ml. So 
we have Si(z) = S2(z)- 

The averaged mass of all black holes is {M,(t)).^n = (M^) + 

/„(M.(f'))aiidf', where (Ml) is the mean mass of seed black 
holes and (M»(f'))aii is the averaged accretion rate of all 
black holes larger than M*. Considering that (M,(f'))aii = 
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6i(t'){M,(t'))qso [see Eq. (9) in Small & Blandford 1992; and 
the first equation in sentence with Eq. (14) in Marconi et al. 
2004], we have (M.(O)aii = {MI) + /,| (fi(f')(M.(f'))qsodf'. On 
the other hand, the duty cycle can be written as S] (f ) = Afact/ Af 
for a single episodic phase within a time interval Af, where 
At = Afact + A/dor; Afact and At dor are the active and the dor- 
mant times, respectively. Then, for the i-th episodic phase 
within Af/, the accreted mass is given by (M,(f/))qso(5i(f/)Af/. 
The mean mass of quasar black holes then reads (M,(f))qso = 
(MS) + ^.(M.(f,'))q,o^i(//)A// = (M^) + /;(M.(f'))qso5i(f')df' 
after multiple episodic phases. We then have (M,(f))aii = 
(M,(f))qso. This relation can be understood more easily if we 
consider multiple episodic growth of black holes. We thus 
have th^ duty cycle of the black holes in a concise form of 

S(z) = di(z) = d2(z). 

The total mass density of the black holes at redshift z con- 
tributed from accretion is then given by 

.dz 



Pacc(z) = 



L.(z) 



V 



-nL,z)dL, (4) 



where L^,(z) is the limit luminosity due to survey sensitivity, 
^(L,z) is the luminosity function, the bolometric luminosity 
is given by LboI = CrLb, with Cb the correction factor relating 
B-band luminosity Lb to LboI, f] the radiative efficiency and c 
the light speed. This density includes all the black holes that 
were brighter than L*(z). The mass density of the active black 
holes at z with L > L* is given by 

/•OO 

Pqso(2)= / HM„z)M,6M„ (5) 

where is the limit black hole mass due to the limit lumi- 
nosity. We apply the Soltan's argument to the present sample at 
any redshift z 



Aot(z) 



roo 



UiM„z)M,dM,=p,,,iz). 



(6) 



Though we do not know the distribution M{M, , z), we know the 
mass density of all black holes from Soltan's argument. Finally, 
we have duty cycle 5(z), 

yOqso(z) 



(7) 



Pacc(z) 

If we know ^{L,z) and $(M,,z), we can easily get the duty 
cycle of quasars at redshift z- We have to stress that Eq. (7) nei- 
ther needs an initial condition nor the assumption of a constant 
Eddington ratio, something which must be specified for in the 
alternative approach solving the continuity equation (e.g. 5 = \ 
at z = 3 in Marconi et al. 2004). 

3. SAMPLE AND BLACK HOLE MASS FUNCTION 

3.1. Spectrum Analysis 

The largest quasar sample, given by Richards et al. (2006) 
from the SDSS DR 3, consists of 15,343 quasars from z = 
to 5, which is complete and homogeneous for apparent mag- 
nitude i= 15-19. We only use those quasars (11,954) with 
z < 2.1 in the present paper. We subtract the continuum and 
iron emission (based on the iron template derived from 1 Zw 
1 spectra, which is kindly provided by R. J. McLure from pri- 
vate communication 2006). We then fit H/3 and Mg II lines. 
For those with z < 0.7, four components are used to model the 
spectra: broad and narrow H/3 plus narrow [O iii] 4959A and 



[O III] 5007A. For others, we use one broad and one narrow 
component to model the Mg ii line. Some objects are removed 
from the sample for one of three reasons; they have: 1) too 
poor spectra to fit due to low signal-noise ratio; 2) only nar- 
row lines (< 2000 km s"') or 3) serious absorption at the Mg 
II line. Finally we have 10,979 quasars available for estimating 
black hole masses with z < 2.1. The removed quasars reduce 
the completeness of the sample given in Table 1 . The last bin 
(zio) is poor because the quahty of quasar spectra is not good 
enough to measure width of Mg II. A future paper will give a 
detailed description of estimating black hole masses and related 
issues (Chen et al. in preparation). 

3.2. Black Hole Mass Function 

We apply the latest version of the empirical relation of re- 
verberation mapping (Kaspi et al. 2005; Vestergaard & Pe- 
terson 2006) to calculate the black hole mass in each quasar. 
For low redshift quasars with z < 0.7, we use the full-width- 
half-maximum of H/? whereas we use Mg ii for those with 
0.7 < z < 2.1 (McLure & Dunlop 2004). The scatter of the 
re-calibrated relation for the black hole masses is less than 0.4 
dex, which is much improved. Figure la shows the mass distri- 
bution of the present sample. We fit the mass distribution of the 
black holes via least square method in a form of three power 
laws 



^(m.) = /o»C 1 + — 1 + — 



m.i 



m,2 



(8) 



where m, is mass of black holes in unit of solar mass. This ex- 
pression has the limits !F{m») oc m^^' for m, <C m,i, Tim,) oc 
^ai-a2 foj. ^ ^ and T{m,) cx m^'""2-"3 for 

m. » OT.2. We obtain /o = (2.70 ± 0.35) x lO'^^, m.i = 
(4.14±0.56) X lO'', m.2 = (2.50±0.22) x 10^ ai = 3.56±0.28, 



ai = 2.87 ± 0.25 and as = 2.71 ± 0.07. A significant break ap- 
pears at M, = 2.5 x IO'Mq and then a steeper mass spectrum 
jF(m,) oc OT,^ '^^ follows the break mass, which is consistent 
with the maximum mass from the SDSS DRl (McLure & Dun- 
lop 2004). 



Table 1 The mass function 



z-bin 




Pi 


Pi 


Ml 


X'j/d.o.f. 


C 




(io-«) 






(lO^Mo) 


(%) 


Zl 


48.1 


-0.60 


1.55 


0.01 


4.40/4 


98 


Z2 


4.74 


0.32 


2.48 


0.16 


7.90/9 


95 


Z3 


5.94 


0.79 


2.58 


0.18 


8.77/10 


93 


Z4 


0.79 


0.18 


2.89 


0.58 


9.89/10 


94 


Z5 


0.44 


0.48 


3.36 


1.25 


8.25/10 


95 


Z6 


0.27 


0.15 


3.21 


1.81 


8.16/9 


96 


Z7 


0.29 


0.42 


2.82 


1.53 


8.78/8 


97 


Z8 


0.33 


0.68 


2.86 


1.55 


8.58/8 


96 


Z9 


0.28 


0.50 


2.54 


1.29 


6.45/7 


89 


ZlO 


0.15 


0.63 


2.30 


1.26 


6.71/7 


67 



The last column is the completeness at each redshift bin in our sample. 
The redshift bin is defined by zj = U~ l/2)Az, where j = 1 , 10. 
The fittings are not included one point in zi, Z4, Z6 and zg -zio-bins, 
which significantly deviate from the trend of double power-laws. 

We get the MF by dividing our sample into 10 redshift bins 
with an interval Az = 0.21 and then into 20 black hole mass 
bins in each redshift bin. Figure lb shows the MF. We find 
that the function can be well fitted by double power laws in the 
following form 



$(M.,z) = $* 



M, 



Pi 



M* 



ft 



(9) 



Duty Cycle of Quasars 
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where /3i, [32, <&* and M* are constants. The peak mass is 
given by M^'' = (/3i//32)^^*^'^'''^m:. We get the four pa- 
rameters from the least x-square method via minimizing = 

X); [^{M, , z) - ^ij] ^ /(jfj, where $y is the MF derived from the 
sample in mass M,.,-bin and redshift z —bin. Since the aver- 
aged error bar in the luminosity function is about A4'/^' ~ 0.1 
(Richai-ds et al. 2006), we take ajj = 0.1$,; for the MF. Table 1 
lists all the parameters for each redshift bin. 



I I 1 1 III 1 — I I I 1 1 III 1 — I I I 1 1 III 1 — I I I 1 1 III 



(a) 




Fig. 1. — a shows the mass distribution of the black holes in our sample. 
The red line is the least square fit. The black hole mass function <I>(M. ,z) at 
each redshift bin is shown in h. The lines represent the best fits given in Tab 1 . 

The MF shows an increasing peak mass toward high red- 
shifts. The mass break at the low mass side is not real. It is 
caused by the survey flux limit. We have to point out that the 
MF is not complete for all black holes, but it is complete for our 
flux-limited criterion. 

4. RESULTS 

The B-band luminosity is converted from the absolute mag- 
nitude Miiz = 2) via a relation of Mb = Mi(z = 2) + 0.804 
(Richards et al. 2006). It has been found that the factor Cb is not 
a function of redshift (Steffen et al. 2006). A more elaborate 
treatment gives Cb = 6-7 for quasars (Marconi et al. 2004). We 
use Cb = 6.5 throughout the redshift range in this paper The lu- 
minosity function is taken from Richards et al. (2006). We take 
from the minimum mass of the black holes in each redshift 
bin to calculate their mass density, pqso(Eq. 5). This corre- 
sponds to the survey limit and is consistent with the accretion 
density pacc (Eq. 4). The radiative efficiency ry > 0.1 is reached 
by Yu & Tremaine (2002), 77 > 0.15 by Elvis et al. (2002) and 
77 = 0.15 by Gammie et al. (2004), based on the Soltan's ar- 
gument, the X-ray background and numerical simulations, re- 
spectively. Numerical calculations indicate that black holes are 
rotating with their maximum spin all the time during their evo- 



lution if accretion is included (Volonteri et al 2005), which is 
supported by studies of SDSS quasars (Wang et al. 2006); then 
77 « 0.3 without significant evolution. We thus calculate the 
duty cycle for two different radiative efficiencies rj = QA and 
0.3. 




Redshift (2) 

Fig. 2. — a, b show the duty cycle of quasars as a function of redshift in our 
sample and the history of the SFR density, respectively. The SFR density is 
taken from Perez-Gonzalez et al. (2005). 

Results are shown in Figure 2. First, quasars have duty cycle 
of 5{z) = 10"^ ~ 1 as shown in Fig 2a, indicating that black holes 
are undergoing active and dormant phases, namely episodic ac- 
tivity. The value of 5{z) ^ 1 at z 2 agrees with the assumption 
of 5 = \ at z = 3 in Marconi et al. (2004). Second, the duty cy- 
cle is rapidly evolving from z « 2 to the local Universe. We 
find (S(z) oc z'^ , where 7 ^ 2.5 until z = 1 .5, above which it tends 
to flatten. Third, as shown in Fig. 2a and 2b, the duty cy- 
cle connects quite naturally with the history of the SFR density 
as a consequence of co-evolution of galaxies and black holes. 
Massive, gas-rich mergers account not only for most of the star 
formation at z « 2-3, but they are probably also responsible 
for triggering major episodes of black hole activity (Di Mat- 
teo et al. 2005). Figures 2a and 2b show that the duty cycle 
follows the star formation history, implying evidence for in- 
trinsic star-formation-triggered black hole activity. The higher 
the SFR density, the higher the triggering frequencies of the 
black hole activity. Last, we define TZ = Afact/Afdor, then have 
7^ = 5{z) / [1 - ^(z)] ■ Episodic activity of the black holes can be 
described by this parameter. When 5{z) 1, black holes have 
TZ ^ 1, namely, the dormant black holes are frequently trig- 
gered at z 2 by star formation. At that time, quasars look like 
long-lived phenomena because of 7?^ 3> 1. 

5. DISCUSSION AND SUMMARY 

The evolution of quasars is jointly controlled by the trigger- 
ing mechanism and accretion. The duty cycle is a key parameter 
to unveil the evolution of quasars. The results of the present pa- 
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per show a very strong cosmological evolution of quasar's duty 
cycle. The triggering history represented by S(z) is quite similar 
to the evolution of cosmic SFR density. This indicates that star 
formation may be the direct mechanism to trigger the activity 
of the black holes. 

The duty cycle can be roughly justified from the galaxy lumi- 
nosity function. According to the luminosity function of galax- 
ies at 1.8 < z < 2.0 (Dahlen et al. 2005), the galaxy number 
density is no ~ 826Gpc"^ for galaxies brighter than /?-band 
magnitude Mr = -24. This corresponds to the number density 
of galaxies with black hole mass larger than lO^M© converted 
from log (Mbh/Mq) = -0.5Mr-3 (McLure & Dunlop 2001). 
Number density of quasars brighter than M,- = -28 (correspond- 
ing to a black hole with mass > lO'^ Mq if it is accreting at the 
Eddington limit) is «q « 175 Gpc"^ based on the quasar lumi- 
nosity function (Richards et al. 2006). We thus estimate a duty 
cycle of ~ 0.18, which is roughly consistent with the present 
results (see Figure 2a). 

The SFR density rises with redshift out to z = 1 -5 and appears 
to be roughly flat between z w 1.5 and z » 3.0. This tendency 



could be explained by strong feedback from activity of black 
holes to their host galaxies (Silk 2005; Di Matteo et al. 2005; 
Croton et al. 2006). With a balance between star formation 
and feedback in z ~ 1 .5 - 3.0, the violent star formation is then 
suppressed. However, the SFR density is going to decrease with 
time due to a shortage of gas and the BH duty cycle follows this 
trend. To further confirm this, future work will focus on the de- 
pendence of the duty cycle on the black hole masses. It could 
show the feedback-dependence on the BH growth itself. Ad- 
ditionally the total accretion time (net Ufetime) of black holes 
will be then obtained. 



The anonymous referee is greatly acknowledged for pointing 
out a mathematical error in the paper The authors are grateful 
to M. Z. Kong, R. Wang and R. J. McLure for help in SDSS data 
analysis, and L. C. Ho, S. N. Zhang, J. X. Wang and X. Y. Xia 
for useful discussions. This research is supported by the Nat- 
ural Science Foundation of China through NSFC-10325313, 
10233030 and 10521001. 



REFERENCES 



Croton, D. J., et al., 2006, MNRAS, 365, 11 

Dahlen, T., et al., 2005, ApJ, 631, 126 

Elvis. M., Risaliti, G., & Zamorani, G., 2002, ApJ, 565, L75 

Gammie. C. R, Sharpiro, S. L., & McKinney, J. C, 2004, ApJ, 602, 312 

Di Matteo, T., Springel, V., & Hernquist, L., 2005, Nature, 433, 604 

Kaspi, S., et al., 2005, ApJ, 629, 674 

Marconi, A., Risaliti, G., Gilli, R., Hunt, L. K., Maiolino, R., & Salvati, M., 

2004, MNRAS, 351, 169 
Martini, P., 2004, in Coevolution of Black Holes and Galaxies, from the 

Carnegie Observatories Centennial Symposia, by Cambridge University 

Press, Ed. L. C. Ho, p. 169 
Martini, P. & Schneider, D., 2003, ApJ, 597, L109 
McLure, R., & Dunlop, J. S., 2001, MNRAS, 327, 199 
McLure, R., & Dunlop, J. S., 2004, MNRAS, 352, 1390 



Perez-Gonzalez. P G., et al. 2005. ApJ, 630, 82 

Rees, M. J., 1984. ARA&A. 22, 471 

Rees, M. J., 1990, Science, 247, 817 

Richards, G. T., et al., 2006, AJ, 131, 2766 

Richstone, D., et al., 1998, Nature, 395A, 14 

Silk, J., 2005, MNRAS, 364, 1337 

Small, T.A., & Blandford, R. D., 1992, MNRAS, 259, 725 

Soltan, A., 1982, MNRAS, 200, 115 

Steffen, A. T., et al. , 2006, AJ, 131, 2826 

Vestergaard, M., 2004, ApJ, 601, 676 

Vestergaard, M., & Peterson, B. M., 2006, ApJ, 641, 689 

Volonteri, M., Madau, P, Quataert, E., Rees, M. J., 2005, ApJ, 620, 69 

Wang, J.-M., Chen, Y.-M., Ho, L. C, & McLure, R. J., 2006, ApJ, 642, Llll 

Yu, Q., & Tremaine, S., 2002, MNRAS, 335, 965 



